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Abstract 

This  paper  presents  a  rigorous  mathematical  formulation  for  modeling  the  upper 
extremity  that  is  capable  of  considering  a  relatively  large  number  of  degrees  of  freedom, 
thus  yielding  a  realistic  model  and  associated  envelope.  Kinematic  models  are  used  to 
determine  the  reach  envelope  in  closed-form  and  to  better  understand  human  motion. 
Joint  ranges  of  motion  are  taken  into  account  by  transforming  unilateral  inequality 
constraints  into  equalities  that  are  included  in  the  formulation.  Methods  from  geometry 
are  implemented  to  analyze  the  motion  and  delineate  barriers  within  the  workspace.  It  is 
observed  that  these  barriers  are  indeed  surfaces  where  the  limb  has  one  or  more  joints  at 
their  limits,  but  also  where  the  hand’s  motion  has  encountered  a  kinematic  singular 
configuration.  Such  a  configuration  is  mathematically  defined  and  is  physically 
associated  with  two  links  being  parallel  at  an  instant  in  time  or  where  two  joints  have 
their  axes  parallel  (e.g.,  a  fully  extended  arm  yields  a  singular  configuration).  Barriers  to 
motion  can  now  be  characterized  in  terms  of  different  human  performance  measures,  thus 
leading  to  a  better  understanding  of  the  path  trajectories  assumed  by  humans  as  they 
execute  tasks. 

Keywords:  Reach  envelope,  ergonomics,  shonlder  complex,  elbow,  wrist. 
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1.  Introduction 

Analytical  methods  in  the  field  of  robotics  have  significantly  contributed  to  the 
evaluation,  programming,  and  design  of  mechanical  robotic  manipulators  and  similar 
devices.  The  motivation  for  investigating  the  motion  of  human  limbs  using  such 
mathematical  methods  is  twofold:  to  achieve  a  better  understanding  of  functionality,  in 
terms  of  evaluative  measures,  and  to  obtain  more  rigorous  methods  for  designing 
ergonomic  workspaces.  The  ultimate  goal  of  this  research,  however,  is  to  develop  a 
method  for  identifying  injured  joints  from  a  path  trajectory.. 

Biomechanical  models  of  the  shoulder  have  been  proposed  by  many  researchers,  such  as 
qualitative  planar  kinematic  models  (Dvir  and  Berme  1978;  Jackson,  et  al.  1977),  but 
were  confined  to  a  single  motion  pattern.  A  complex  model  of  the  shoulder  was 
proposed  by  Hogfors,  et  al.  (1991),  in  which  it  was  treated  as  a  three-rigid  body  twelve 
degree-of-freedom  (DOF)  system.  From  a  biomechanical  point  of  view,  the  human  arm 
mechanism,  particularly  the  shoulder  joints,  are  probably  the  most  complex  of  the  human 
body.  The  upper  extremity  is  typically  modeled  as  a  series  of  segmental  links  connected 
in  a  special  arrangement  of  one  DOF  revolute  or  prismatic  joint  (Wood  et  al.  1989, 
Benati  et  al.  1980  and  Engin  et  al.  1990).  These  links  have  typically  been  limited  to  a 
small  number  of  DOFs  because  of  the  difficulty  in  addressing  larger  numbers.  We  will 
follow  a  similar  modeling  method,  but  will  not  limit  the  number  of  DOFs.  Instead,  we 
will  develop  and  demonstrate  a  mathematical  method  for  analyzing  human  motion, 
visualizing  the  resulting  workspace  to  better  understand  the  barriers  therein.  Most 
importantly,  we  will  account  for  ranges  of  motion  in  determining  and  visualizing  the 
envelope. 

In  an  earlier  investigation  (Lenarcic  et  al.  1994),  a  simple  model  of  a  human  arm  only 
considers  the  shoulder  complex  and  the  elbow  joint  where  the  two  translational  motions 
were  replaced  with  two  rotational  DOFs  to  simplify  the  analysis.  The  European  Esprit 
Project  CHARM  (Maurel  1998)  developed  a  comprehensive  human  animation  resource 
model  allowing  the  dynamic  simulation  of  complex  musculoskeletal  systems,  including 
finite  element  deformation  of  soft-tissues  and  muscular  contraction.  A  biomechanical 
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model  of  the  human  arm  and  shoulder  was  designed  including  properties  for  bones, 
joints,  and  muscle  lines  of  action.  However,  reach  envelopes  were  not  determined. 
Indeed,  many  researchers  (including  commercial  human  modeling  software  developers) 
have  adapted  a  simplified  approach  of  moving  the  arm  through  its  ranges  (arm  fully 
extended)  to  determine  the  outer  most  boundary.  While  this  method  yields  only  one 
barrier,  it  is  approximate  in  nature,  does  not  determine  all  barriers,  and  does  not  yield 
closed  form  equations  of  the  boundary  surfaces. 

This  paper  extends  our  previous  work  (Abdel-Malek  &  Yang  2001)  by  presenting  a 
practical  model  of  the  human  arm  that  includes  the  resulting  motion  of  the  shoulder 
complex.  The  first  section  of  this  paper  reviews  the  modeling  approach  and  provides  a 
brief  anatomy  of  the  upper  extremity.  In  the  second  section,  we  present  the  kinematic 
model  and  the  theory  of  swept  volumes  that  we  will  use  to  determine  the  reach  envelope. 
Based  on  this  approach,  we  use  two  steps  to  analyze  the  human  arm  workspace.  First,  we 
will  directly  address  the  workspace  of  the  7DOFs  without  the  translational  motion  of  the 
scapulathorasic  joint.  Second,  we  will  subject  the  resulting  geometry  to  a  sweep  motion 
due  to  the  shoulder  translational  joints.  We  are  then  able  to  visualize  the  reach  envelope, 
characterize  boundary  surfaces  in  closed  form,  and  describe  barriers  therein  that  are  due 
to  joint  limits. 


2.  Background 

In  order  to  develop  a  realistic  kinematic  model  of  the  upper  extremity,  it  is  important  to 
understand  human  anatomy.  The  upper  limb  is  composed  of  three  chained  mechanisms: 
the  shoulder  girdle,  the  elbow  and  the  wrist.  If  we  consider  bones  in  pairs,  seven  joints 
may  be  distinguished:  the  stemo-clavicular  joint,  which  articulates  the  clavicle  by  its 
proximal  end  onto  the  sternum,  the  acromio-clavicular  joint,  which  articulates  the  scapula 
by  its  acromion  onto  the  distal  end  of  the  clavicle,  the  scapulo-thoracic  joint,  which 
allows  the  scapula  to  glide  on  the  thorax,  the  gleno-humeral  joint,  which  allows  the 
humeral  head  to  rotate  in  the  glenoid  fossa  of  the  scapula,  the  ulno-humeral  and  the 
humero-radial  joints,  which  articulate  both  ulna  and  radius  on  the  distal  end  of  the 


4 


humerus,  and  finally,  the  ulno-radial  joint  where  both  distal  ends  of  ulna  and  radius  join 
together. 


The  shoulder  girdle  is  perhaps  the  most  diffieult  to  model  because  of  the  extension  of 
some  muscles  over  more  than  two  segmental  links  and  joints.  The  two  main  bones  of  the 
shoulder  are  the  humerus  and  the  scapula  (shoulder  blade).  The  end  of  the  scapula,  called 
the  glenoid,  meets  the  head  of  the  humerus  to  form  a  glenohumeral  cavity  that  acts  as  a 
flexible  ball-and-socket  joint  (Figure  1). 


Figure  1 .  Anatomy  of  the  shoulder  complex 

Shoulder  movements  are  usually  referred  to  as  ventral/dorsal,  cranial/caudal  and  axial 
rotations  for  the  stemo-clavicular  (3  DOFs),  as  abduction/adduction,  flexion/extension 
and  axial  rotation  for  the  gleno-humeral  joint  (3  DOFs),  as  elevation/depression, 
protraction/retraction,  tipping  forward/backward  and  medial/lateral  rotations  for  the 
scapulo-thoracic  joint  (5  DOFs),  and  as  flexion/extension  and  pronation/supination 
movements  for  the  forearm  joints  (2  DOFs).  In  our  modeling,  we  shall  use  SDOFs  for 
the  shoulder  complex,  where  we  have  accounted  for  3  rotational  and  two  translational 
motions  (Figure  2). 
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Figure  2.  Modeling  of  the  shoulder  complex  as  three  revolute  and  two  sliding  DOFs 

This  model  allows  for  consideration  of  the  coupling  between  some  of  the  joints,  as  is  the 
case  in  the  shoulder  where  muscles  extend  over  more  than  one  segment.  When  muscles 
are  used  to  lift  the  arm  in  a  rotational  motion  a  translational  motion  of  the  shoulder 
unwittingly  occurs. 

The  elbow  is  a  hinge  joint  made  up  of  the  humerus,  ulna  and  radius.  The  unique 
positioning  and  interaction  of  the  bones  in  the  joint  allows  for  a  small  amount  of  rotation 
and  hinge  action.  Theprimary  stability  of  the  elbow  is  provided  by  the  ulnar  collateral 
ligament,  on  the  medial  (inner)  side  of  the  elbow.  It  is  safe  to  model  the  elbow  as  a  IDOF 
revolute  joint. 

The  hand  is  composed  of  many  small  bones  called  carpals,  metacarpals,  and  phalanges. 
The  two  bones  of  the  lower  arm  —  the  radius  and  the  ulna  —  meet  at  the  hand  to  form  the 
wrist.  We  will  model  the  wrist  as  three  revolute  joints  intersecting  at  one  point,  whose 
action  yields  a  spherical  wrist  (Pieper  1968).  The  complete  9  DOF  model  of  the  upper 
extremity  is  shown  in  Figure  3. 
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Figure  3.  Modeling  of  the  upper  extremity  as  a  9  degree  of  freedom  kinematic  chain 

Given  the  respective  material  properties,  bones  may  be  regarded  as  rigid  bodies  in 
contrast  to  soft  tissues  with  respect  to  the  relevant  physiological  ranges  of  motion  and 
force  handling.  This  allows  for  a  rigid  body  segmental  representation  of  the  kinematics 
among  the  skeletal  subsystem  from  the  soft  tissues  by  converting  their  relations  with  the 
bones  into  external  actions.  Therefore,  the  kinematic  model  can  only  be  analyzed  by 
considering  the  skeletal  components. 

3.  Mathematical  Formulation 

Although  the  anatomy  of  the  human  upper  extremity  and  joints  are  very  complex  as 
evidenced  by  the  debate  in  the  literature  on  the  correct  method  for  modeling  joint  motion 
(Engin  1984;  Choi  1993;  Engin  1989a,b;  Hogfors  1991;  Raikova  1992;  Helm  1994, 
Maurel  1998),  we  will  employ  a  kinematic  pair  (or  combination  thereof)  used  in  the  field 
of  robotics.  In  fact,  all  anatomical  joints  can  be  modeled  using  a  combination  of  basic 
kinematic  pairs.  Consider  the  two  segmental  links  connected  by  a  joint  in  Figure  4.  We 
will  use  a  coordinate  frame  for  each  degree  of  freedom  in  the  system.  Because  the 
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shoulder  in  our  model  has  5DOFs,  for  example,  it  will  have  five  variables  assoeiated  with 
it  in  whieh  eaeh  variable  is  denoted  by  q. 


Figure  4.  Definition  of  a  kinematic  pair  (e.g.,  a  revolute  joint) 

We  define  q*  H[^i  ...  *(55iR"  as  the  vector  of  n-generalized  coordinates 

characterizing  the  motion  of  a  limb,  where  qi  (/  =1,  2,  ...  ,  n)  represents  a  DOF.  For 
instance,  the  elbow  joint  will  be  represented  by  one  DOF,  for  example  qs  (Figure  4), 
while  the  wrist  joint  will  be  represented  by  three  DOF,  namely  qy,  qs  and  qg.  The  vector 
function  generated  by  a  point  of  interest  written  as  a  multiplication  of  rotation  matrices 
and  position  vectors  is  expressed  by 

(1) 

1=1  7=1 


0(q*)  = 


x(q’) 

T(q) 

2(q) 


where  both  '  'p,.  and  ^  are  defined  using  the  Denavit-Hartenberg  (DH) 
representation  method  (Denavit  and  Hartenberg  1955),  such  that 


COS0. 

-COS  a 

sin  0. 

sin  « .  sin  0. 

-R,  = 

sin  0^ 

cos  a. 

cos  0. 

-  sin  cos  6. 

(2) 

0 

sin 

cos  a^ 

and 

(i- 

II 

o 

o 

a-  sin0j 

(3) 
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where  0.  is  the  joint  angle  from  the  x._j  axis  to  the  x.  axis,  J,.  is  the  shortest  distance 
between  x._j  and  x.  axes,  a,  is  the  offset  distance  between  z,.  and  z,._j  axes,  and  a,  is 
the  offset  angle  from  z,_j  and  z,  axes.  The  generalized  variable  is  =  0.  for  a 
revolute  joint  and  q.  =  d.  for  a  prismatic  joint. 


The  vector  function  0(q)  characterizes  the  set  of  all  points  inside  and  on  the  boundary 
of  the  reach  envelope  generated  by  an  anatomical  landmark,  and  typically  selected  as  a 
fingertip.  The  objective  is  to  visualize  this  vector  function  consisting  of  many  parameters 
and  to  better  understand  the  motion  governed  by  ®(q*) .  At  a  specified  position  in  space 
given  by  P(x^,y^,z^) ,  Eq.  1  can  be  written  as  a  constraint  function. 


n(q) 


x(fl)-Xp 

T(q)-T, 

z(q)-Zp 


=  0 


(4) 


Ranges  of  Motion 

Ranges  of  motion  are  imposed  in  terms  of  inequality  constraints  in  the  form  of 
q-<q,<q-,  where  /  =  1,  ...  n,  are  transformed  into  equalities  by  introducing  a  new  set 
of  generalized  coordinates  k  =  [/I,  ...  such  that 


where  the  new  variable  is  inherently  constrained  by  the  sine  function  and  does  not 
change  the  dimensionality  of  the  problem.  In  order  to  include  the  effect  of  the  ranges  of 
motion,  it  is  proposed  to  augment  the  constraint  equation  Q(q*)  with  the  parameterized 
ranges  of  motion,  such  that 


H(q)  = 


x(q’)-x^ 

yi^)-yp 

q.  -  a.  -  b.  sin  X- 


=  0  1  =  1,  ...  n 


(6) 
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where  q  = 


^ 

q  A 


■^T 


T  i 

is  the  vector  of  all  generalized  coordinates.  Note  that  although 

n-new  variables  (1.)  have  been  added,  equations  have  also  been  added  to  the 
constraint  vector  function  without  loosing  the  dimensionality  of  the  problem. 


The  Jacobian  (after  the  German  mathematician  Karl  Gustav  Jacob  Jacobi)  of  the 
constraint  function  H(q)  at  a  point  q®  is  the  (3  +  n)  x  2n  matrix 

H^=^H/^q  (7) 

where  the  subscript  denotes  a  derivative.  With  the  modified  formulation  that  includes 
ranges  of  motion,  the  Jacobian  is  expanded  to 

O.  I  O’ 

H,  =  (8) 

where  q^  =  5q*/5X, ,  =  53>/ dq’ ,  0  is  a  (3  x  n)  zero  matrix,  I  is  the  identity  matrix, 

and 


0)  .  = 

q 


A 


z 


^1 


^2 

y,.  - 


^2 


% 


-((^r-^f)/2)cos/i, 


0 

0 

0 


0 

-((^2^  -^2  )/2)cos/i2 
0 
0 


0 

0 

0 


-((q::-q:)l2)cos\ 


(9) 


(10) 


In  earlier  work  by  the  authors  aimed  at  determining  difficulties  in  the  control  of  robot 
manipulators  (Abdel-Malek  and  Yeh  1997a;  1997b;  Abdel-Malek,  et  al.  1997;  and 
Abdel-Malek,  et  al.  1999),  it  was  shown  that  impediments  to  motion  (halting  of  a 
trajectory)  arise  inside  the  workspace  when  the  Jacobian  becomes  singular  (Abdel-Malek 
and  Yeh  2000).  Because  the  Jacobian  is  not  square,  rank  deficiency  criteria  were 
developed.  Before  addressing  these  criteria,  it  is  important  to  show  why  the  singularity 
of  the  Jacobian  has  a  direct  effect  on  the  control.  The  differentiation  of  Eq.  (1)  with 
respect  to  time  yields  the  velocity  of  the  fingertip  O  as 

®  =  (11) 
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where  q*  is  the  vector  of  joint  velocities.  Given  a  specified  path  trajectory  (i.e.,  ® ),  the 
calculation  of  q*  (i.e.,  joint  velocities)  requires  computing  an  inverse  of  the  Jacobian  . 

For  a  singular  Jacobian,  it  is  not  possible  to  compute  the  required  velocities.  These  cases 
are  typically  associated  with  a  kinematic  configuration  of  the  upper  extremity  that  does 
not  admit  motion  in  a  particular  direction,  but  requires  a  change  in  the  arm’s  posture  in 
order  to  execute  the  path.  If  the  Jacobian  was  square,  then  the  determinant  of  will 

yield  the  postures  in  space  where  singular  behavior  occurs.  We  will  use  this  concept  to 
explore  the  surrounding  workspace. 

Jacobian  Analysis 

Because  the  Jacobian  is  not  square,  we  define  these  singular  sets  as  a  subset  of  the 
workspace  in  which  the  Jacobian  of  the  constraint  function  of  Eq.  (7)  is  row  rank 
deficient  (Abdel-Malek,  et  al.  1999);  i.e.,  the  barriers  are  defined  by  and 
characterized  by 

^c{RankHq(q)  <k,  for  some  q  with  H(q)  =  0} 

where  k  is  at  least  (3  +  i7-l).  Imposition  of  the  rank  deficiency  condition  can  be 
implemented  using  a  variety  of  methods,  but  perhaps  the  most  computationally  efficient 
one  is  the  repeated  elimination  of  square  sub-Jacobians,  until  a  number  of  non-linear 
equations  are  determined.  For  example,  consider  a  7  DOF  model  of  the  arm,  where  the 

and 

3x7 


Therefore,  three  sets  of 


l7x7 

singularities  can  be  identified: 


Jacobian  is  10x14,  where  is  in  the  following  form. 


where  q^  is  in  this  form: 


(12) 
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(a).  The  largest  square  submatrix  of  is  3x3  -  --  --  -- 

JsxV 


n\  7! 

-------  There  are  a  total  of - = - =  35  submatrices  (also 

3!(«-3)!  3!4! 

■^35  J3x7 

called  sub-Jacobians)..  The  determinants  of  these  matrices  are  |jj|=-  - 

I/35I  =  -  -  -  .  If  all  the  determinants  are  kept  at  zero,  then  the  group  of  equations  can 
then  be  numerically  solved  to  identify  the  singular  sets. 

(b).  When  q.  reaches  its  limit,  eliminate  the  /  th  column  in  and  it  will  be  a  3x6 

matrix.  Therefore  there  are  20  sub-Jacobians  and  20  equations  can  be  solved  together  to 
obtain  the  third  set.  When  q.  and  q.  reach  their  limits  repeat  the  same  procedure  to  find 

the  singular  set.  When  q^ ,  qj  and  qj^  reach  their  limits  repeat  the  same  procedure  to  find 

the  singular  set  until  the  remaining  Jacobian  of  is  a  square  matrix.  These  sets  are 

denoted  by  a,  b,  and  c,  where  each  set  comprises  a  number  of  constant  joint  values  (from 
the  total  number  of  DOF)  and  two  variables.  Upon  substituting  these  sets  into  Eq.  (1),  we 
obtain  an  equation  of  a  surface  parameterized  in  terms  of  two  variables  (i.e.,  a  2DOF 
surface  in  3D  space),  such  that: 

f('>(u<'))  =  ®(u<'),qO  (13) 

where  — >■  q*  fl  and  f is  the  vector  function  describing  the  new  parametric 

surface  characterizing  a  barrier  to  motion. 
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(c).  Since  is  a  square  matrix  the  singular  sets  can  be  obtained  by  solving  the  equation 


* 


=  0 .  The  singular  sets  are  exact  when  the  joints  reach  their 


limits. 


In  order  to  separate  the  analysis  of  rotational  joints  from  that  of  translational  joints,  we 
break  the  formulation  down  into  two  distinct  steps.  First,  we  address  the  all-revolute  7 
DOF  model  of  the  human  upper  extremity  (including  glenohumeral,  elbow  and  wrist),  but 
exclude  the  two  shoulder  translational  joints.  We  will  then  perform  a  sweep  of  the 
resulting  envelope  to  visualize  the  complete  reach  envelope. 


Motion  on  the  Barriers  of  Motion 

To  better  understand  when  the  hand  may  or  may  not  cross  barriers  under  given 
conditions,  we  explore  the  barrier’s  kinematic  properties.  We  propose  a  criterion  that  is 
based  on  normal  acceleration  at  a  point  on  a  barrier,  such  that  crossability  is  achieved  if 
the  barrier  admits  a  normal  acceleration  in  one  direction  or  another.  A  point  on  a  barrier 
admits  motion  normal  to  the  surface  in  either  direction  depending  on  the  difference  in 
acceleration  components  (defined  by  the  indicator//),  such  that 


7  = 


Po 


(14) 


where  v,  is  the  tangential  velocity,  is  normal  acceleration,  and  Ij is  the  normal 
curvature  of  the  barrier  with  respect  to  the  tangent  direction  of  v,  {p^  is  the  radius  of 
curvature).  The  need  for  formulating  the  problem  in  terms  of  velocities  and  accelerations 
will  become  apparent,  as  the  resulting  expression  for  the  indicator  rj  will  be  independent 
of  acceleration  values,  but  will  be  a  quadratic  form  that  has  definiteness  properties.  A 
point  on  a  singular  surface  will  have  no  acceleration  if  the  quantity  rj  computes  to  null. 
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For  a  singular  parametric  entity  e  (where  u  is  a  vector  representing  the 

remaining  joint  variables — those  not  constant),  and  in  the  field  of  differential  geometry, 
the  First  Fundamental  Form  (Farin  1993)  is  denoted  by  I^,  where 

v]  =1^;  qj\  ,  and  is  defined  as 

(15) 

where  f„  =  da..  The  Second  Fundamental  Form  is  defined  as 

II  da  (16) 

P  L  Juu 

or  expanded  to 

II^  =  Wf^Ju^  +mxju  dv  +  Wf„dv^  (17) 

where  N  is  the  vector  normal  to  the  singular  surface  and  =  ^ijdud^ .  The  Normal 
Curvature  of  a  parametric  surface  at  a  configuration  q^,  in  the  direction  of  dujdv , 
can  then  be  defined  as  the  ratio  (Farin  1993) 

^:,  =  — =‘C  (18) 

Po  Ip 

In  order  to  determine  the  kinematics  quantities,  we  define  the  Time-Modified  First  and 
Second  Fundamental  Forms  as 

i;  =  uXf„u  (19) 

ir^u^fN^fl  u  (20) 

P  L  Juu 

such  that  the  normal  curvature  can  still  be  defined  as 


p  I  r 

I  o  p  p 


(21) 


For  a  singular  surface  f^'^(u^‘^),  the  derivative  using  the  chain  rule  is  f„u.  Similarly,  for 
the  overall  description  of  the  workspace  G(q),  the  derivative  is  G^q.  Therefore,  at  an 

instant  of  time,  the  tangential  velocity  in  terms  of  f  or  G  at  any  point  on  the  barrier  is 

V,  =  f>  =  G,q  (22) 
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If  joint  limits  are  considered,  then  the  derivatives  can  be  written  as  u  =  u^s  and  q  =  q^s . 
The  squared  norm  of  the  velocity  is 

\  \2  T  •  TrTr  • 

K|  =V,V,  =u  f„f„u 

which  is  equal  to  the  Time-Modified  First  Fundamental  Form  of  Eq.  19.  Therefore, 
Ij,  can  be  written  as 


(23) 


Substituting  1/p^  into  rj  yields 


(24) 


I  f 


(25) 


This  expression  can  be  written  in  quadratic  form  whereby  the  matrix  of  the  quadratic 
form  need  only  be  evaluated. 


4.  Reach  Envelope  of  the  Upper  Extremity  Excluding  Transiation 

Consider  a  model  of  the  upper  extremity  comprising  7  rotational  joints  (7DOFs)  as 
shown  in  Figure  5.  Note  that  this  model  does  not  include  the  two  translational  joints  in 
the  shoulder,  which  will  be  addressed  in  the  subsequent  section.  Note  too  that  the  arm  is 
fully  extended,  as  it  was  in  the  initial  configuration.  We  will  consider  the  following 
ranges  of  motion  from  this  initial  configuration: -90°  <  ^3  <  90° ,  -1 10°  <  ^4  <  120° , 

-90°  <  ^3  <90°,  -150°<^g<0°,  -60°  <  ^3  <60°,  -20°<^8<20°,  and 

-90°  <^9  <90°. 
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Tracked  Point 


Figure  5.  Kinematic  model  of  the  7DOF  all  revolute  upper  extremity 


According  to  the  DH  method,  the  XYZ  coordinates  of  the  tracked  fingertip  point  in  terms 
of  the  seven  variables  are 

®(q)[l]  =  -lOc^c^s^  +  20s^s^  + 1 5cg  (-c^c^s^  +  )  - 1  + 1 0(Cg  (c^  (Cg  {-c^c^s^ 

+  >^3^5  )  -  )  +  (^5*^3  +  C3^4‘^5  )^7  )  +  "  {c^C.S^  +  )5g  )58 

®(q)[2]  =  -IQc^S^S^  -  2OC3 ^5  +  1 5Cg  (-Cji'j ^4  -  C3 ^5  )  - 1 504^35'^  +  1 0(Cg  {Cj  (Cg  (-C5 535'4 

-  )  +  (-C3C5  +  ^354^5  )s,  )  +  (-C4Cg53  -  (-C55354  -  C355  )5g 

0(q)[3]  =  2OC4C5  + 1 5c4C5Cg  - 1 5s^s^  + 1 0(Cg  (Cj  (c^c^c^  -  ^4^6 )  -  04^553 ) 

+  (-Cg54  -C4C5^g>8) 


Singularity  Sets 

n\  7! 

Since  ^>  .  is  a  3x7  matrix,  there  are  - = - =  35  equations  to  be 

^  3!(«-3)!  3!4! 

simultaneously  solved,  which  represent  the  determinants  of  all  square  3<l>3  sub- 

Jacobians.  There  are  three  sets  of  solutions: 

(a).  As  a  result  of  solving  the  equations  generated  by  the  sub-Jacobians,  the  rank 
deficiency  of  yields  Sj  =  {q^  =  0,  ^7  =  0,  =  0} .  This  set,  when  substituted  into  Eq.  1 


yields  a  surface  parameterized  into  three  generalized  coordinates. 


(26) 
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(b).  The  second  set  is  identified  when  one,  two  or  three  joints  reach  their  limits  (in  this 
example,  ®(q)  only  includes  6  variables,  When  one  joint  reaches  an  upper  or 

lower  limit,  for  example  q^=2^/3 ,  we  substitute  q^=2^/3  into  0(q)  and  compute 
(the  remaining  Jacobian  excluding  the  q^  column).  Since  is  a  3x6  matrix, 
6! 

there  are  =  20  equations  to  be  simultaneously  solved,  whose  solution  is 
a  =  {^5  =  0,  q-i  =  0,  =  0} .  Therefore,  a  singular  set  is  identified  as 


'^71 

=  {^4  =  =0,^7  =0,^g  =0}.  Substituting  only  one  of  the  constraint  limits  and 


applying  the  rank-deficiency  condition  yields  the  singular  sets  listed  in  Appendix  A.  This 
process  continues  for  all  variables.  When  two  variables  reach  their  constraint  limits,  such 
as  q^=-7t  I2,q^=27t  13 ,  substituting  these  limits  into  ®(q)  and  calculating  (again. 


excluding  both  associated  columns)  yields  4  equations.  When  simultaneously  solved  we 
obtain  b  =  {^7  =  0,  =  0} .  Therefore  a  singular  set  is  identified  as 
1*2  =  {^3  =  =  2;7r/3,^7  =  0,^g  =  0} .  Repeating  the  process  also  yields  the 


singular  sets  listed  in  Appendix  A.  When  three  variables  reach  their  constraint  limits  such 


as  ^3  =  -y,^4  =  ~  ’  ®^b®tituting  these  limits  into  0(q)  and  calculating 


(again,  excluding  three  associated  columns)  yields  1  equation, 
c  =  {^8  =  0} .  Therefore  a  singular  set 
1  Xk 


n 


*^1  {^3  2  ’^4 


n 


-,q.  = - ,^g  =0}. 

18  '  2  ' 


When  solved,  we  obtain 
is  identified  as 


(c).  The  third  set  is  when  four  of  the  variables  reach  their  constraint  limits  such  as 

-;^/2l  r-ll;T/18l  \-7rl2\  \-57rl6\ 

W2  )•  2;r/3  f  =U/2  f  =1  0  f 

singular  surfaces  in  Appendix  A.  In  order  to  perform  this  symbolic  manipulation,  we 
have  developed  a  computer  code  using  Mathematica®. 
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The  determinants  of  sub-Jacobians  for  ®(q)  are  as  follows: 


l”'.l  =  KI  =  -  =  KI  =  o 

\Jn\  =  500((2  +  cos^7  cos^g)sin^g  +cos^g  sin^g)(-cos^g  sin^4  sin^^  sin^^  +cos^4((2 

+  cos q-j  cos ) sin q^  +  cos q^  sin q^)  +  cos q^  sin q^ (4  +  cos q^ (2  +  cos q.j  cos q^)- sin q^  sin q^ )) 


IJ29I  =  |•>^3o|  =  250  cos  (2  cos  ^5  cos^2  +sin^5  sin  ^2  )((2  +  cos  ^2  cos^g)sin^g  +cos^g  sin^g) 

I/jjI  =  -250 cos ^g (cos ^g  sin^j  sin^g  sin^^  +  cos ^5 (2 cos ^7  cos^g  sin^g  +sin^g  +2cos^g  sin^g)) 

IJ32I  =  250  cos  ^g  sin  ^7  (4  +  2  cos  ^7  cos^g  +  cos  ^g  (2  + cos  ^7  cos  ^g)- sin  ^g  sin^g) 

I/33I  =  125sin^g(2(l  +  cosgg)(cos^7)^  cos^g  +(-3  +  cos(2^7))cos^g  +  2  cos  ^7  (4 +  2  cos  ^g  -sin^g  sin^g)) 
IJ34I  =  -250(2  + cos  ^gXcos^g)^  sin  q-, 

1/35!  =  -125sin(2^g) 

Enforce  all  above  to  be  zero  and  solve  the  equations  simultaneously  and  one  obtains  the 
rank  deficiency  of  which  is  Sj  =  {^g  =  0,  ^7  =  0,  ^g  =  0} . 


When  q.  reaches  its  limit,  such  as  ^3  =-— ,  we  substitute  q^  =-nl7.  into  0(q)  and 

compute  O^.  (the  remaining  Jacobian  excluding  the  q^  column).  The  determinants  of  the 
sub-Jacobians  are 
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|c/i|  I^/jI  ■■■  l-^iol  0 

I^jjI  =  -500((2  +  cos^7  cos ^g) sin +cos^g  sin^g)(-cos^g  sin ^5  sin ^7 
+  cos ^5 (4  +  cos (2  +  cos cos q^)- sin ^g  sin q^ )) 

=  -250(1 +  2  cos  ^g)  cos  sin^5(-cos^g  sin^j  sin  ^7  +  cos  ^5  (4  +  cos  ^g  (2  +  cos  ^7  cos^g) 

-sin^g  sin^g)) 

IJjjI  =  -250(2  cos  ^g  sin^g  +  (l  +  2cos^g)cos^7  sin^g)(-cos^g  sin^g  sin  ^7  +  cos  ^3  (4  + cos  ^g 
(2  +  cos  ^7  cos  ^g )  -  sin  q^  sin  q^ )) 

1714!  =  250 cos  ^g (2 cos  ^5  cos  ^7  +sin^5  sin^7)((2  +  cos^7  cos  ^g)  sin  ^g  +cos^g  sin^g) 

I^jjI  =  250(cos^7  sin  ^5  -2  cos  ^5  sin  ^7 )  sin  ^g  ((2  + cos  ^7  cos  ^g)  sin  ^g  +cos^g  sin^g) 

l^jgl  =  -125(2(cos^g)^  sin^g  sin^g  sin^7  +  cos  ^5  (4  cos  ^7  (cos  ^g)^  sin^g  +(l  +  2cos^g)sin(2^g))) 

1717!  =  250  cos  ^g  sin  ^7  (4 +  2  cos  ^7  cos^g  +  cos  ^g  (2  + cos  ^7  cos  ^g)- sin  ^g  sin^g) 

l^igl  =  125sin^g(2(l  +  cos^g)(cos^7)^ cos^g  +  (-3  +  cos(2^7 )) cos ^g  +2cos^7(4  +  2cos^g  -sin^g  sin^g)) 

|jj9|  =  -125sin(2^g) 

\J20\-  -250(2  + cos ^gXcos^g)^  sin ^7 


If  we  enforce  the  rule  that  all  the  determinants  will  be  zero  simultaneously  and  solve 

TT  TC 

these  equations,  we  obtain  one  singular  set,  ag  ={^3  =  —,q^  = - ,q^  =0,^8  =0}  .  One 


can  then  obtain  other  a^  by  repeating  this  procedure.  When  q^  and  q .  reach  their  limits. 


TT  1  XtT 

such  as  ^3  = - and  q^  = - ,  we  substitute  them  into  ®  (q)  and  compute  ®  .  (the 

2  18  ** 


remaining  Jacobian  excluding  the  q^  and  q^  columns).  The  determinants  of  the  sub- 
Jacobians  are 

kihk2h---=Kho 

\j.j  \  =  250 cos  ^8  sin  ^7  (4  +  2  cos  ^7  cos^g  +  cos  ^g  (2  + cos  ^7  cos  ^g)- sin  ^g  sin^g) 

|jg|  =  125sin^8(2(l  +  cos^g)(cos^7)^  cos^g  +(-3  +  cos(2^7))cos^8  +  2  cos  ^7  (4 +  2  cos  ^g  -sin^g  sin^g)) 
[/gl  =  -250(2  + cos  ^g)(cos^g)^  sin  ^7 
j/jgl  =  -125sin(2^g) 
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Once  again,  we  enforce  all  of  the  determinants  to  be  zero  simultaneously  and  solve  these 

71  1  XtI 

equations  to  obtain  one  singular  set  bj={^3= - ,^4= - ,^3=0,^g=0}.  By 

2  18 

repeating  this  procedure  other  are  obtained.  When  ,  qj  and  q^.  reach  their  limits, 
71  1  ItT  71 

such  as  ^3  = - ,  q^= - and  q^= - ,  we  substitute  them  into  0(q)  and  compute 

2  18  2 

(the  remaining  Jacobian  excluding  the  q^ ,  q^  and  q^  columns).  The  determinants  of 

sub-Jacobians  are 

kil=k2i=k3i=o 
\Ja\  =  -125sin(2^g) 


We  enforce  all  the  determinants  to  be  zero  simultaneously  and  solve  these  equations  to 

71  1  XtI  71 

obtain  one  singular  set  Cj  =  {^3  = - ,^4= - - repeat  this 

2  18  2 


procedure  to  obtain  other  . 


When  q*  is  singular,  then  q* 


1 1 5.7r^  cos  /I3  cos  /I4  cos  cos  cos  cos  X^  cos  /I5 
419904 


=  0. 


We  can  obtain  X^  and  plug  it  into  Eq.  5;  that  means  q^  reaches  its  limit.  While  0(q) 
doesn’t  include  q^  when  any  four  of  the  variables  reach  their  limits,  it  will  be  singular 

set.  It  is  now  possible  to  visualize  each  surface  by  substituting  all  singular  sets  into  Eq. 
26,  as  illustrated  in  Figure  6. 
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^169  "^184 


^201  "X2I6 


Figure  6.  The  different  singular  sets  of  surfaees 
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In  fact,  combining  all  singular  surfaces  yields  the  reach  envelope  of  the  7  DOF  model  as 
shown  in  Figure  7.  We  have  identified  the  outer  boundary  to  the  reach  envelope  in 
closed  form  (also  shown  in  Figure  8). 


Figure  7.  The  eross  seetion  view  of  workspace  for  human  arm  model  (7  DOF)  and  eross 

seetions 

at  z  =  0  and  z  =  20 
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Figure  8.  Envelope  of  the  7  DOF  Model  for  Human  Model  and  the  cross  sections 

at  z  =  0  and  z  =  20 


5.  Sweeping  the  7-DOF  in  two  translational  directions 

The  effect  of  the  scapula-thorasic  translational  motion  is  accounted  for  using  a  swept 
volume  analysis  (Blackmore,  et  al.  1999).  We  now  include  two  translational  joints  to  the 
shoulder  complex,  where  the  kinematic  model  of  a  9  DOF  upper  extremity  is  shown  in 
Figure  9. 
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1 


Figure  9.  Kinematic  model  for  9  DOF  Human  Arm 


There  are  a  number  of  surface  patches  that  constitute  the  outer  boundary.  These  surfaces 
will  be  denoted  by  F,(q),  where  m  is  the  number  of  surface  patches,  and 

q  =  has  no  more  than  two  variables  for  each  surface  patch.  A  sweep  operation 

in  the  direction  of  ^(^1,^2)  =  [^'i  ^2  translational  directions)  with  no 

rotation  is  defined  by  the  sweep  equation 

=R  0,(q)  +  y{q„q2) 

Swept  volume  Swept  entity  Translational  path 

where  the  rotation  matrix  R  is  I ,  an  identity  matrix  since  no  changes  occur  in  its 
orientation,  =  \q^,q^,u,v^  now  contains  four  variables  for  each  surface  (Figure  10), 

where  u  and  v  are  the  parametric  variables  of  the  surface  patch  on  the  boundary  of  the  7 
DOF  reach  envelope. 


(27) 
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Figure  10.  Sweeping  of  the  7DOF  outer  boundary  in  two  translational  directions 

In  order  to  visualize  this  vector  function  r(p)  ,  we  again  make  use  of  the  rank  deficiency 
condition,  where  the  Jacobian  5r/0p  is  a  3x4  matrix.  The  four  determinants  of  the 
sub- Jacobians  of  5r/5p  are  set  to  zero  and  singular  sets  are  determined.  As  an  example, 

35cos(w)cos(v) 

consider  the  sweep  of  the  surface  patch  defined  by  ®(m,v)=  35cos(M)sin(v)  and 

35sin(v) 

characterized  by  the  equation 
-i-35cos(m)cos(v) 

r(p)=  ^2 +35cos(M)sin(v)  ,  where  0<u<l7r,  -6<v<nll,  -1.5<^j<1.5  and 
35sin(v) 

-1.5  <  ^2  ^  1-5 . 
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The  Jacobian  matrix  is  Fp 


1  0  -35cos(v)sin(M)  -35cos(M)sin(v) 

0  1  -35sin(i/)sin(v)  35cos(m)cos(v) 

0  0  0  35cos(v) 


.  Singular  sets 


have  been  determined  with  the  following  values:  g^={u  =  Q,v  =  7r  12} , 
g^  =  {u  =  n,q^=-\.5] ,  g^=  {u  =  Q,q^=\.5}  ,  g^  =  {u  =  'in  l2,q^=-\.5] , 

g^={u=7i:  I2,q^=\.5} ,  and  g^  =  {v  =  Q,q2=-\.5} .  These  sets  are  automatically 
substituted  into  r(p)  to  visualize  the  boundary.  The  process  is  repeated  for  all  surface 
patches,  and  the  resulting  reach  envelope  is  shown  in  Figure  1 1 . 
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20  0  -20  -40 


Figure  1 1 .  The  final  envelope  of  9  DOF  human  arm  model 
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Barriers  to  Motion 


To  better  understand  the  relation  between  barriers  within  the  reach  envelope  and  the 
physical  meaning  of  a  singular  configuration,  we  examine  a  singular  set  characterized  by 

a  =  =0,^2  =0,^3  =0,^g  = - =  -7r/3,q^  =  -7r/9,qg  =0)5  shown  in  Figure  12 

6 


TT  K 

with  the  configurations  {q^= - ,q^=0  ,q^=  —  ).  In  this  case,  there  is  no  translational 

18  12 


Stt 

motion  of  the  scapula  thorasic  motion;  however,  the  elbow  is  at  the  limit  q^  = - ,  and 

6 

the  wrist  joint  has  reached  a  combination  of  limits  (^7  =  0,^g  =  -;r/9,^5  =0).  Figure  13 
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shows  the  configurations  of  three  singular  sets  ^5  =  -— ,  ^5  =  0 ,  and  ^5  “  • 
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Figure  12.  Configurations  of  three  singular  sets 
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71  71 

F  igure  1 3 .  Three  eonfigurations  combined  together 


The  Jacobian  analysis  has  yielded  seven  constant  values  of  the  joint  variables,  but  has  left 
two  variables  unspecified  (^4,^5).  Substituting  the  set  of  constant  variables  into  Eq.  27 
yields  a  new  equation  parameterized  in  terms  of  the  two  remaining  variables  as 


r(^4>^5) 


I  (8  -  3^/3  +  2  cos(^)  -  2  Vs  sin(^))  sin(^5 ) 

I  (cos(^4 )  cos(^5  )(-8  -  3  Vs  +  2  cos(^)  -  2  Vs  sin(^))  -  (3  +  2  V3  cos(^)  +  2  sin(^))  sin(^4 ) 
I  (cos(9,  )(3  +  2  V3  cos(^)  +  2  sin(^))  +  co^{q,  )(8  -  3  V3  +  2  cos(^)  -  2  V3  sin(^))  sin(^4 ) 


Note  that  the  above  expression  represents  the  sweep  of  the  fingertip  in  space  with  all 
joints  but  q^,q^  in  a  locked  position.  We  plot  this  2DOF  surface  in  3D  space,  as  shown 
in  Figure  14. 
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Figure  14.  A  barrier  due  to  singular  set 
Stt 

2^  =  {q^  =0,^2  =0,^3  =0,^g  = — —,q^  =  -nll>,q^  =-nl9,q^  =0} 

6 

The  physical  significance  of  this  surface  can  be  observed  by  examining  the  inverse 
kinematic  solution  on  any  point  on  the  surface.  This  particular  configuration 
characterizing  the  barrier  delineates  important  properties  of  the  reach  envelope. 

To  explain  this  physical  behavior  typically  associated  with  Jacobian  and  limit 
singularities,  consider  the  motion  of  one's  arm  aimed  at  reaching  (with  the  fingertip)  a 
point  on  the  barrier  surface  (Figure  7),  where  many  joints  have  reached  their  limits  (or 
have  reached  a  Jacobian)  singularity,  and  some  joints  are  free  to  move.  For  a  given 
configuration  with  the  fingertip  at  this  point,  it  is  not  possible  to  admit  motion  normal  to 
the  barrier  in  at  least  one  direction.  However,  the  barrier  is  indeed  in  the  reachable 
workspace;  therefore,  any  point  to  the  left  and  right  of  this  point  are  reachable.  This 
means  that  a  different  configuration  (other  than  the  configuration  that  yields  the  barrier), 
may  allow  the  fingertip  to  admit  motion  across  the  barrier.  Such  a  result  has  also  been 
observed  with  mechanical  manipulators  (especially  in  the  field  of  welding),  where 
continuous  motion  of  the  end-effector  is  necessary.  If  the  weld  seam  is  positioned  across 
a  barrier,  the  manipulator  may  not  smoothly  execute  the  motion,  but  may  halt,  change 
inverse  solutions,  and  then  continue  the  motion.  We  believe  that  such  kinematically- 
driven  behavior  of  the  upper  extremities  is  important  in  identifying  disabilities  in  issuing 
motor  commands,  and  could  play  an  important  role  in  identifying  such  disabilities. 


6.  Conclusions 

A  9  degree  of  freedom  biomechanical  model  of  the  upper  extremity  has  been  developed, 
and  was  used  to  generate  the  reach  envelope.  We  have  presented  a  rigorous  formulation 
for  modeling,  analysis,  and  visualization  of  the  envelope  and  barriers  where  one  or  more 
joints  have  reached  their  limits  or  where  a  Jacobian  singularity  is  identified. 
Fundamental  to  this  formulation  is  the  underlying  concept  of  accounting  for  joint  ranges 
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of  motion  in  terms  of  inequality  constraints  imposed  on  the  motion.  It  was  shown  that 
rank  deficiency  conditions  of  the  Jacobian  of  the  model  yield  surface  patches  that  belong 
to  the  envelope.  It  was  also  shown  that  a  subset  of  these  surfaces  fall  on  the  boundary. 
Translational  joints  of  the  shoulder  are  accounted  for  using  a  swept  volume  approach, 
which  yields  closed  form  equations  of  the  boundary.  Furthermore,  barriers  within  the 
reach  envelope  are  important  in  explaining  and  verifying  hand  trajectories.  Barriers  that 
are  a  result  of  kinematic  singularities  present  obstacles  to  motion  for  the  upper  extremity 
when  encountered  during  a  path,  requiring  that  the  motion  be  halted  and  switched  to 
another  posture.  We  believe  that  the  identification  of  these  barriers  will  aid  in  better 
understanding  motor  commands  issued  by  the  central  nervous  system,  as  well  as  perhaps 
providing  insight  into  path  following  for  people  with  disabilities. 

This  study  is  the  first  step  towards  a  more  rigorous  investigation  of  the  reachable 
workspace;  a  forthcoming  extension  of  this  work  is  to  include  upper  extremity  orientation 
capabilities  to  better  understand  human  dexterity.  Also  underway  is  a  study  to  relate 
barriers  that  appear  within  the  workspace  to  neurological  commands  and  motor  skills. 
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9.  Appendix  A 

27r  1  Itt 

^1  ={^'4  =0’^8  =0}’  ^2  =  {^4  =  --—,^5  =0,^7  =0,^3  =0}, 

J  i  o 

71  1  XtI  71  1  XtI  71 

^3  “  {^3  “  “  ^  ?  ^4  “  9  ^5  “  “  ^  ’  ^8  “  ’  ^4  “  ’  ^5  “  ^  ?  ^7  “  ^9  ^8  “  9 

^5  “  {^3  “  ^  9  ^5  “  “  “  9  ^7  “  O9  ^8  “  9  ^6  “  {^3  “  ^  9  ^5  “  “  9  ^7  “  O9  ^8  “  9 

—  —  =0,^8  =0},  83  ={^4  =  —,^5  =-y, ^2  =0,^3  =0}, 

=-^^1  =09^8  =0}^  ={^4  =0^^8  =0}^ 

^,,={q,=-^,q,=0,q,=0,q^=0}  ,  =  {q,  =  ^,q,  =  0,q,  =  0,q,  =  0}  , 

^13  “  {^4  “  “^9^6  “  ^9^7  —~~->q^  “  0}  9  ^14  —  {^4  —  “9^6  “  ^9  ^7  “  ““9^8  “  ’ 

^15  “  {^5  “  ^9  ^6  “  ^7  “  “^9  ^8  “  ’  ^16  “  {^4  ~  “  ^9  ^7  “  “9  ^8  “  ’ 

^17  “  {^4  “  ^  9  ^6  “  O9  ^7  ~  ~->q^  “  0}  9  ^18  “  {^5  “  ^9  ^6  “  O9  ^7  ~  ~->q^  “  0}  9 

TT  1  ItT  TT  StT 

bi  =  {73  =-^^^4  =  -- ^^77  =0,78  =0}>b2  ={73  =--,^4  =--,72  =0,^3  =0}, 

2  io  2  3 

TT  1  ItT  TT  StT 

^^3  “  {^3  “^9^4  “  .  ^  9  ^7  “  O9  ^8  “  ^}  9  1^4  “  {^3  ~  ~Z  ~  ^  9  ^7  “  O9  ^8  “  ^}  9 

2  io  2  3 

*>5  ={73  =-y.75  =-y.77  =0.78  =0}.  bg  ={73  =-y.75  =-y.76  =  0,78  =  0}, 

^7  “  {^3  “  ^9  ^5  —  ~~->qq  “  O9  ^8  “  9  bg  —  {^3  —  ~~->q^  “  “~9  ^6  “  O9  ^8  “  9 

^9  “  {^3  ~  “  O9  ^8  “  ’  ^10  “  {^3  ~  ~  “  ^9  ^8  “  ’ 

bn  ’^5  =^9^7  =0^  ^8  =0}’  bi2  ={^3  =^9^6  =0^^8  =0}^ 
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^13  “  {^3  “  ~  ’  ^14  “  {^3  “  ~~^->^5  ~  ~~^->^6  ~^‘>^1  “  ? 

^15  “  {^3  ~  ~~^->^5  ~~^->^6  “  0?  ^7  “0}  ’  ^16  “  {^3  ~~^->^6  “  0?  ^7  “  0?  ^8  “  ’ 

^17  “  {^3  “  ^5  ~  “  0?  ^7  “  ^}  5  ^18  “  {^3  “  ^5  “  ^6  “  0?  ^7  “  ’ 

^19  “  ~~^->^5  ~~^->^6  ~  ^’^7  “  ’  ^20  “  {^3  ~~^->^5  “  “  ^’^7  “  ^}  ’ 

2  2  6  2  2  6 

^21  “  {^3  “  ~  ~  ^’^7  “  ^}  ’  ^22  “  {^3  “  ~  ~  ^’^7  “  ^}  ’ 

2  2  6  2  2  6 

b23  ={^3  =0,^7  =-y,^8  =0}^  *>24  =  {^3  =  ^ ,q,  =  0,q,  = -^ =  0}  , 

*>25  ={^3  =0^^7  =y.^8  =0}^  *>26 

1  Itt  tt  tt 

•>27  ={^4  =  -- r^^76  =0>^7  =-^’^8  =^^’'^28  =3,^6  =0,^7  =--,^8  =0}, 

1 0  J  J 

1  ItT  TT  27r  TT 

*>29  “  {^4  “  J  g  ’  ^6  “  0?  ^7  “  ^8  “  ^}  ’  *>30  “  {^4  “  ^  ’  ^6  ~^->^l  “  ^8  “  ^}  ’ 

*>31  “  {^5  “  ~~^->^6  ~  ^9^8  “  ^}  ’  *>32  “  {^4  “  ~~->^5  ~  ~~^->^6  ~^->^l  ~  0} 

*>33  “  {^4  “  ^9  ^5  “  “^9  ^6  “  O9  ^7  “  0}  9  *>34  “  {^4  “  O9  ^5  “  “~9  ^6  “  ^9  ^7  “0} 

*>35  “  {^5  “  ^  9  ^6  “  O9  ^7  “  O9  ^8  “  ^}  9  *>36  “  {^4  “  “  ^  9  ^5  “  ~  9  ^6  “  ^9  ^7  “0} 

*>37  ={^4  =^>^5  =^9^6  =0>^7  =0}>  *>38  ={^4  =0>^5  =  ^  9  ^6  =  0>  ^8  =  0} 

b39  ={^4  =-^,^5  =-^,^6  =-^>^7  =0}^  *>40  ={^4  =^>^5  =-^>^6  =-^>^7  =0} 

2  2  6  2  2  6 

641  ={^4  =-^,^5  =^,q,=-^,q,=0},  b,,={q^=^,q,  =^,q^=-^,q^  =0} 

2  2  6  2  2  6 

*>43  ={75  =-y.76  =0.77  =-y.78  =0}>  644  =  {^5  =0,^7  =-y,^8  =0} 

*>45  ={75  =-y.  76  =0,77  =y,78  =0},  64^  ={^5  =  0,^7  = -y,^8  =  0} 


TT  1  Itt  tt  tt  1  Itt  tt 
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TT  27r  TT  TT  1  ItT  TT 
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TT  27r  TT  TT  1  ItT  TT 
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^6  ^  -arccos( - - - ^ ^ - )} 
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2(- 1 6  +  8  cos  ^0+3  cos  q^  -  a/-48  sin  q^  +104  cos  q.  sin  q^  -35  cos  q^  sin  q^  +  4 sin  q.^ ) 

q^  =  -arccos( - - r - r - )} 

16  +  8cosgg +cosgg  +4singg 

,  TC  WtT  7Z 

C32  ={^3  =“Y’^4  =-  — .^7  =-> 

2(-16  +  8coso„  +3cosff/  +  J-48sinfl/  +104cosff„  sino/  -35coso/  sino/  +4sino/) 

=  -arccos( - - - - r - )} 

16  +  8cos^g +cos^g  +4sin^g 

,  71  3.71  71 

^33  “  {^3  “  ^  ?  ^4  “  ’  ^7  ~  ~  ^  ’ 

2(- 16  +  8  COS  ^g +  3  cos  q^  -  J-4S  sin  ^g^  + 1 04 cos  q.  sin  q^  -35  cos  q^  sin  q^  +  4 sin  q^ ) 

q^  =-arccos( - - - r - r - )} 

16  +  8cos^g +cos^g  +4sin^g 

,  TT  371  7t 

=y.^4  =  — ’^7  =--’ 

2(-16  +  8coso„  +3cosff„^  +  J-48sinfl„^  +104cosff„  -35coso7  sinfl„^  +4sino/) 

=  -arccos( - - - - r - )} 

16  +  8cos^g +cos^g  +4sin^g 

2(-16  +  8cos^g  +3cos^g^  - ^-48 sin q^  +104 cos q^  sin q^  -  35 cos q^  sin q^  +  4 sin q^ ) 


,  7t  WtC  7t 

^35  “  {^3  “  ^  ?  ^4  “  “  ’  ^7  “  ^  ’ 


q^  =  -arccos( 


16  +  8  cos  ^g  +cos^g^  +4sin^g 


)} 


,71  WtI  7C 

C36={73=y.^4=-  — ^^7=-. 

2(-16  +  8coso„  +3cosff„^  +  J-48sinfl„^  +104cosff„  sino/  -35coso7  sinfl„^  +4sino/) 

q^  =-arccos( - - r - r - )} 

16  +  8cos^g +cos^g  +4sin^g 

TT  371  7C 

C37  =  {^3  =  “  ~  ?  ^4  =  -7^  ?  ^7  “  ’ 

2(- 16  +  8 cos  ^g +  3  cos  q^  -  J-4^ sin  ^g^  + 1 04 cos  q.  sin  q^  -35  cos  q^  sin  q^  +  4 sin  ^g"^ ) 

q^  =  -arccos( - - r - r - )} 

16  +  8cos^g +cos^g  +4sin^g 

,  TT  37T  7T 

C38=fe=-Y’74=— .^7=-’ 

2(-16  +  8coso„  +3cos7„^  +  J-48sino„^  +104cos7„  smo„^  -35coso7  sinOj,^  +4sino/) 

q^  =  -arccos( - - - - r - )} 

16  +  8cos^g +cos^g  +4sin^g 

2(-16  +  8cos^g  +3cos^g^  -  ^-48  sin  q^  +104 cos  q^  sin  q^  -  35  cos  q^  sin  q^  +  4  sin  q^ ) 


.  TT  37T  7T 

C39  =  {^3  =  ~  ?  ^4  =  -7^  ’  ^7  “ 


q^  =  -arccos( 


16  + 8  cos  ^g  +cos^g^  +4  sin  q^ 


)} 
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,  71  ^71  71 

^40  ~  {^3  “  ^  ?  ^4  “  ’  ^7  “  ^  ’ 


=  -  arccos( 


2(-16  +  8cos^g  +3cos^g^  +y-48sin  + 104 cos ^g  sin  -35  cos ^g^  sin^g^  +4sin^g'^)^ 

16 +  8  cos  ^g  +cos^g^  +4sin^g^ 


C41  ={^3  =-y’^6  =0’^7  =-J’^8  =0}’  C42  ={^3  =y’^6  =-y’^8  =0} 

C43  ={^3  =-^’^6  =-^’^7  =-^’^8  =^^’^44  ={^3  =-^’^6  =^’^7  =^’^8  =0^ 

Z  O  J  Z  J 

y4 

C45  —  {^3  —  ~  9  ^6  “  ^  ’  ^7  “  ~  ^  5  ^8  “  ^}  9  ^46  “  {^3  ~~Z->^6  ^  ?  ^8  “ 

2  6  3  2  3 

c,,={q,=-^,q,=-^,q,  =^^q^=0},  c,,  =  {q,  =^^q^=-^^q^  =^^q^=0} 

2  6  3  2  6  3 

^49  “  {^3  “0?^?  ““  ^}  9  ^50  “  {^3  “  ?  ^5  “  ^  ?  ^6  “  0?  ^7  “  “ 

.  ;r  ;r  5;r  ;r.  tt  tt  ^  tt  . 

^51  “  {^3  ~~7->^5  ~  ~  ^’^7  “  ~~2^  ’  ^52  “  {^3  ~~7->^5  ~  “^?^6  ~  “} 

4  2  6  3  4  2  3 

C53  =^^^5  =^^^6  =-^^^7=-^}^  ^54=  {^3  =^^^5  =^^^6  =^^^7  =  ^} 

4  2  6  3  4  2  3 

.  TT  71  StI  71.  .  71  71  StI  71. 

^55  “  {^3  “  ^9^5  “  “^9^6  “  ^’^7  “  9  ^56  “  {^3  “  ^9^5  “  ~9^6  “  ^’^7  “ 

4263  4263 

.71  7Z  „  TT-,  .71  7Z  „  TT. 

C57  =fe  =^’75  =--.^6  =0>^8  =-^}7  C58  ={^3  =-,^5  =-,^6  =0,^8  =-^} 

.  ;r  ;r  5;r  ;r.  .  tz  tz  ^  ti  ^ 

C59  ={73  =^’75  =-y»76  =-— .78  =--}7  Cgo  =  {^3  =-.75  =-y.^6  =0,^8  =-^} 

.  7Z  StZ  7Z .  7Z7Z  7Z 

^61  “  {^3  “  9  ^5  —  ~  9  ^6  “  9  ^8  “  “  '^}  9  ^62  “  {%  ~~->^5  “  ^  9  ^6  “  ^9  ^8  “ 

.  StZ  7Z  .  .  7Z  7Z  StZ  7Z  . 

^63  ={73  =^.75  =-^.76  =  — ^.78  =^}»  C64  ={73  =^.75  =“.76  =  —^.78  =  ^} 


("  A  ^  ^3  ("  57r  TT  TT , 

ri  =  {76  =0,77  =-T’^8  =“7r^’  ^2  =  {76  =-— .77  =-^’78  =--} 

3  9  6  3  9 

73  ={76  =0,77  =y,78  =-^}.  74  =  {76  =  0,77  = -y.78  =^} 

.  5;r  ;r  ;r.  .  5;^  ;r 

75  =  {76  =  -^’77  =  ^.78  =  -^-l .  76  =  {76  =  -^.77  =  - ^.78  =  -8-} 

6  3  9  6  3  9 

77  ={76  =  0,77  =-.78  =-}.  78  ={76  =—r^7i=-,qi=-} 

3  9  6  3  9 

r-;T/2l  r-ll;7r/18l  [-ttH 

79. ..724  Sixteen  combinations  q^=\  [’^4=1  f’  "  I  jtH 


and^g  = 


-Stt  !  6 


37 


^25  •  •  •  Yw  Sixteen  combinations  = 


-nil) 


741 .. .  Xjg  Sixteen  combinations  ^3  = 


-n  19 


Y^i-.. Y-,2  Sixteen  combinations  ^3  = 


-nil 


Y73---Yss  Sixteen  combinations  = 


-n  19 


Yw-Ym  Sixteen  combinations  q^  = 


-n  19 


Ym---Ym  Sixteen  combinations  q^  = 


-nil 


Ym---  Ym  Sixteen  combinations  q^  = 


\-5n ! 6\  \-nl9 

^6=  .  h^8= 


Yni---Y\$2  Sixteen  combinations  ^3  = 


-n  19 


Y\si---Ym  Sixteen  combinations  q^  = 


-nil 


-nil 


2W3  •‘'=“U/2 


-nil 


-nil 


W2  2W3  •‘'=“U/2 


-nil 


-nil 


-1  1;t/18 


C)  Qa=-\  r » =  i 

nil  "  Inll  "  0 


-5n !  6 


-nil 


-1  1;t/18 


C)  Qa=-\  C)  Qf.=\ 

nil  "  Inll  '  0 


-5n !  6 


-nil 


-ll;r/18 


C)  Qa=-\  r 3  ^7  =  i 

nil  "  Inll  ’  nil 


-nil 


-nil 


-nil 


;r/2  ’"‘n  0 


-5n !  6 


-nil 


r »  ^3  =  i 

nil  nil 


-nil 


-nil 


-nil 


r  3  ^3  =  i  r  3  ^7  =  i 
nil  '  ;r/2  ’  nil 


-nil 


-ll;r/18 


-nil 


Inll  '  nil  "  0 


-5nl  6 
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/i69  •  •  •  Tm  Sixteen  combinations 


7i85  . . .  ^200  Sixteen  combinations  q^ 


7201- ••7216  Sixteen  combinations  q^ 


f-llir/lSl  \-7rl2\  [Stt  1 6^ 

1  2;r/3  J’ =  U/2  j’ =1  0  j 


J-lliir/lSl  {-7rl2\  {-ttI?) 

I  Itt I2>  |;ir/2j’^’  |;ir/3 


r-;ir/2l  [Stt  !  6\  [-7113 

U/2r«‘=i  0 


f-lla-/18l 

7217  •••7232  Sixteen  combinations  q^=\  K  ^6  “ 

[  l7rl3  J 


-71 13 

71 13 
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